High-field recovery of the undistorted triangular lattice in the frustrated metamagnet 

CuFe0 2 . 
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Pulsed field magnetization experiments extend the typical metamagnetic staircase of CuFeCb up 
to 58 T to reveal an additional first order phase transition at high field for both the parallel and 
perpendicular field configuration. Virtually complete isotropic behavior is retrieved only above this 
transition, indicating the high- field recovery of the undistorted triangular lattice. A consistent phe- 
nomenological rationalization for the field dependence and metamagnetism crossover of the system 
is provided, demonstrating the importance of both spin-phonon coupling and a small field-dependent 
easy-axis anisotropy in accurately describing the magnetization process of CuFeCb- 

PACS numbers: 75.30.Kz, 75.10.Hk, 75.80.+q, 75.30.Gw 
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Metamagnetism typically refers to any material that, 
upon variation in the externally applied magnetic field, 
exhibits an abrupt change in magnetization. In general, 
the phase diagrams of materials undergoing field-induced 
magnetic transitions can be rationalized according to 
the degree of magnetic anisotropy in the materials^. In 
highly anisotropic systems, spins are effectively restricted 
to align (anti-)parallel to the magnetic easy-axis and 
magnetic transitions typically involve discontinuous spin 
reversals, leading to first-order type metamagnetic tran- 
sitions. As for isotropic (weakly anisotropic) systems 
this directional restriction is relieved (strongly reduced), 
transitions in such materials often reflect the onset of a 
continuous, second order type reorientation of the local 
spins. Another source of exotic magnetic transitions is 
geometrical magnetic frustration, which occurs when a 
specific lattice geometry prevents the simultaneous min- 
imization of all magnetic exchange interactions, thus in- 
troducing a high spin degeneracy 2 -. The simultaneous 
occurrence of both these phenomena and the interplay 
between them leads to intricate, diverse and rich physics, 
yielding many captivating magnetic phases ranging from 
spin liquids and ices to multiferroic spiral phases 3 -! 4 ^ 5 -^. 

Here the focus is on the delafossite semiconductor 
CuFe02, an archetype triangular lattice antiferromagnet, 
in which the Fe 3+ ions stack in hexagonal layers along 
the c-axis (Fig. [2(a)). I n spite of the expected Heisen- 
berg nature of the Fe 3+ spins (3d 5 , S = 5/2, L = 0), 
CuFe02 does not order in the noncollinear 120° spin con- 
figuration at low temperature. Instead, after undergo- 
ing successive phase transitions at Tni ~ 14 K and Tn2 
w 11 K, lowering the symmetry from hexagonal (R3m) 
to monoclinic 7,8,9 , the system adopts a collinear, two- 
up two-down order, with moments aligned (anti) parallel 
along the c-axis (Fig. [2(b)) 10 - The collinear ground 
state is supposedly stabilized through the strong spin- 
lattice coupling in CuFeO^^ 1 ^, which induces a struc- 



tural distortion through the 'spin Jahn- Teller' effec t 12 ! 13 . 
Alternatively, this scalene triangle distortion has been 
argued to induce an easy axis anisotropy, which was 
also used to account for the the observed Ising-like 
magnetis m 14 ! 15 ! 16 . An intriguing behavior arises when 
CuFe02 is subjected to an external magnetic field B be- 
low Tn2- With B\\c, the spin system has been found 
to successively assume a proper helical ordered ferroelec- 
tric phase, a collinear three- up two-down ordered phase, 
a phase with a magnetization plateau at one-third of 
the saturation value and a phase with steadily increasing 
magnetization 5 - ! 15 ! 16 ! 17 ! 18 . Particular attention has gone 
to the ferroelectric helical ordered phase 5 - ! 19 ! 20 ! 21 , which 
has since also been stabilized in zero field through Al 3+ 
or Ga 3+ substitutio n 22 ! 23 ! 24 . With increasing B±c, the 
magnetization has been found to first increase steadily, 
then halt at a one-third plateau before resuming a (quasi- 
) linear increase at higher field a 15 ! 16 . This article presents 
pulsed field magnetization experiments, which extend the 
metamagnetic staircase of CuFe02 up to fields exceed- 
ing 58 T— . An additional high field first order phase 
transition is observed for both parallel and perpendicular 
configurations, above which virtually complete isotropic 
behavior is retrieved, indicating the recovery of the undis- 
torted triangular structure. Moreover, a consistent phe- 
nomenological interpretation is provided, combining all 
magnetic terms deemed to be important in CuFe02- 

As depicted in Fig. [2j the magnetization process for 
both B\\c and B_Lc shows a cascade of phase transitions, 
in excellent agreement with literature 5 - ! 15 ! 16 ! 17 ! 18 ! 26 . As 
B I' increases, the spin system successively rearranges to 
the helical ordered phase at b\ x - 7.2 T (Fig. [2(c)) and 

the collinear three- up two-down phase at b\ 2 — 13.0 T 
(Fig. [2(d)). Recent synchrotron x-ray-diffraction studies 
up to 40 T revealed the strong correlation between the 
spin Jahn- Teller lattice distortion and the magnetization 
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FIG. 1: (Color online) (a) Schematic crystal structure of 
CuFe0 2 (R3m, a=6=3.03A, c=17.17A). To avoid confusion, 
crystal directions are referred to using the hexagonal descrip- 
tion throughout the article, (b-e) Spin structures and lattice 
symmetries in various field-induced phases of CuFeCb (-S||c): 
collinear four-sublattice (4SL) phase ((b)), ferroelectric heli- 
cal phase ((c)), collinear 5SL phase ((d)) and collinear 3SL 
phase ((e)). 




5 10 15 20 25 30 35 40 45 50 55 60 
B(T) 

FIG. 2: (Color online) Pulsed field magnetization process for 
both B\\c (red/light gray) and B±c (blue/dark grey) at 1.5 
K. Inset: Zoom-in on high field region. 



order (Fig. D](e)) 16 i 18 , which is consistent with an ex- 
pected nonzero D and pulsed field nuclear forward scat- 
tering experiments^. At b\± ~ 32 A T, the system un- 
dergoes a second order transition, above which M starts 
growing continuously, indicating gradual canting of the 
spins away from collinearity. By extending the experi- 
mental field of view to higher fields, the persistence of 
this (quasi-) linear increase up to b\§ ~ 53.3 T could be 
determined, where an additional metamagnetic transi- 
tion is identified. Above this transition, up to 58.3 T, M 
grows steadily once more, with slightly different slope. 
In short, below b\ a magnetic transitions are of first or- 
der, exhibiting significant hysteresis and large magneti- 
zation steps. From b\ a on (second order transition), M 
mostly increases continuously with i?", and magnetiza- 
tion plateaus are absent. In terms of metamagnetism, 
this corresponds to a crossover from a highly anisotropic 
regime (with abrupt spin flips) to a weakly anisotropic 
regime (continuous spin reorientation), which is in line 
with the notion of progressive symmetry increase and 
thus magnetic anisotropy reduction upon increasing B. 

For B-Lc, the magnetization process is quite different. 
Starting from the zero-field collinear two-up two-down 
phase, M shows a steady increase up to B^ ~ 24.8 
T, where a first order transition brings the system in 
a plateau phase at one-third of saturation, implying a 
three sublattice (3SL) structure. Above B^ 2 — 30.0 T, 
the system exhibits a steady increase of M after under- 
going a second order phase transition, indicating a con- 
tinuous spin reorientation. The data in Fig. [2] show this 
behavior persists up to B^ ~ 51.6 T where the system 
undergoes an additional first order transition, similar to 
that at B\ b for B\\c. Contrary to previous claims^, the 
behavior clearly remains anisotropic up to these tran- 
sitions. Though left unaddressed, a corresponding fea- 
ture can also be observed around 52 T in the (dM/dB) 
vs. B data previously recorded in a single turn coil mea- 
surement up to 100 T (8 K) 26 . At fields above both 
these transitions, M shows virtually isotropic behavior, 
growing (quasi-) linearly up to 58.3 T. This absence of 
anisotropy suggests a full symmetry recovery and thus re- 
trieval of the undistorted triangular lattice at these fields. 

To elucidate the nature of the high field spin structures, 
we introduce a simple classical spin model for a single tri- 
angular sheet, which includes the primary terms in the 
spin Hamiltonian: 



process in CuFeC^; as M increases, the extent of the 
distortion decreases accordingl y 16 ! 18 . Since the induced 
magnetic anisotropy (D) is directly coupled to the distor- 
tion, one may assume it also diminishes correspondingly 
with B] exhibiting steps at first order transitions and 
continuously decreasing in (quasi-) linear phases. More- 
over, at b\ 2 , the symmetry of the distortion is increased, 
yielding a lattice of isosceles triangles 9 . At B^ 3 ~ 19.7 
T, the magnetization jumps to a plateau phase at one- 
third of saturation, signaling a collinear two-up one-down 



H = -gfj,B -^Si + Yl Ji i Si ' S i 

-^bJ l3 (S l -S J ) 2 -D(B)Y,S^ (1) 

where is the exchange coupling between sites i and j, 
b is an effective (nearest-neighbor only) biquadratic inter- 
action originating from the spin-lattice coupling (bond- 
phonon model)^^ and D(B) is the anisotropy constant 
(> for an easy-axis along z). The first (Zeeman) and 
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last (anisotropy) terms sum over all sites i in the mag- 
netic unit cell, while the summation in the exchange and 
biquadratic terms includes all interactions within that 
unit cell. The high field magnetization process is qual- 
itatively similar for both field configurations; first the 
spin system exhibits a one-third magnetization plateau, 
implying a three sublattice (3SL) structure, after which 
M starts increasing steadily, indicating a continuous re- 
orientation of the 3SL spins. To capture this magnetic 
behavior, we thus study the spin Hamiltonian of the 3SL 
structure on a single sheet. As the lattice distortion per- 
sists up to at least 40 T for both field configurations^, 
a nonzero D may be expected up to these fields. Having 
three inequivalent spins, there are three unique first- (J), 
second- (jt) and third- neighbor (J*) couplings per spin 
in one magnetic unit cell (Fig. [3(a)). The Hamiltonian 
becomes (using g = 2, Si = eiS (unit vector e, classical 
spins) and pij = • ej): 

H = -2fi B SB • * + CS 2 (p 12 +Pi3 +P23) + 9J f 5 2 

i 

- GS\p\ 2 +P 2 13 +P 2 23 ) - D{B)S* £ e?„ (2) 

where the exchange constant C = (3 J + 3 J*) and the 
spin-lattice constant G = 3b J. The anisotropy D(B) 
is approximated to be 'an^-proportional' to M(B); as 
M approaches saturation, D vanishes accordingly. Note 
that the inclusion of further neighbor interactions has 
only a trivial effect; the third neighbor interactions 
merely add to C, while the second neighbor interactions 
only shift the total energy as a whole. 

We determined the spin directions e\ , e 2 and es corre- 




FIG. 3: (Color online) (a) Exchange interactions in the 3SL 
structure on the 2D isosceles triangle lattice. Minimum en- 
ergy solutions of eq. [2] for a 3SL structure with finite D are 
depicted for B\\c ((b)) and B±c ((c)). 

sponding to the minimum energy per unit cell at a given 
field B by performing numerical minimization of eq. El 
In examining the corresponding simulated magnetization 
curves, one finds that having both the anisotropy con- 
stant D(B) and the spin-lattice constant G nonzero is a 
prerequisite for good qualitative agreement with exper- 
iment; with D = the magnetization process becomes 
isotropic, while a nonzero G is required to stabilize a 
magnetization plateau for B_Lc. Matching the simulated 
magnetization curves with experiment yields estimates 
for the exchange, anisotropy and spin-lattice constants 
of C = 1.32 meV, L>(25T) = 0.021 meV (plateau) and 



G = 0.0074 meV, respectively, the resulting curves are 
plotted in Fig. [H Taking only first-neighbor interactions, 
we can estimate JS 2 as - 2.76 meV (32.0 K) and DS 
(at 25T) as ~ 0.052 meV (0.6 K), in line with previous 
estimatesJJ-14»l&. For GS A we estimate - 0.29 meV (3.4 
K), yielding a dimensionless biquadratic coupling b of ~ 
0.0056 (compared to ~ 0.008 using estimates from ref. 
11). With these parameters the simulations are in strik- 
ing agreement with experiment. The spin-lattice inter- 
action G (non-directional) stabilizes the one-third mag- 
netization plateau in both configurations, while the di- 
rectional anisotropy interaction D widens the plateau for 
B\\c and narrows it for £?±c, leading to the difference in 
plateau- widths and above-plateau increase of M. More- 
over, a nonzero G also induces a positive d 2 M/dB 2 in 
the latter, as is observed in experiment. 

According to the minimum energy solution the three 
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FIG. 4: (Color online) Simulated magnetization process (eq. 
E} in CuFe0 2 for both B\\c (thick, red) and B±c (thin, blue). 
M is plotted in the both the field direction (left axis, M") 
and the plane perpendicular to B (right axis, M ). Thick 
dark (thin light) grey lines depict offset experimental data for 
B\\c (B-Lc). Upper inset: field-dependence of z-components 
of individual spins for B\\c. Lower inset: estimated magnetic 
energy gain w.r.t. isotropic spin structure (see text). 

spins evolve as follows for B\\c: below i?J 4 , the system 
is in the collinear 3SL state, with two spins (Si and S 2 ) 
parallel to , and one (S3) antiparallel. As depicted in 
the upper inset of Fig. [H from b\ a on, the 'down' spin 
starts continuously tilting from 'down' to 'up', thereby 
increasing M . In optimizing the overall magnetic energy, 
the two 'up' spins respond by first moving slightly away 
from the c-axis in the opposite direction, before gradu- 
ally returning after the 'down' spin has passed the basal 
plane. Due to the finite D and G, the two 'up' spins re- 
main collinear (Fig. E9(b)) and the system also acquires 
a small in-plane magnetization, which quickly grows and 
slowly decreases with above b\ a (Fig. [4]). Although 
the process is qualitatively similar for B_Lc (Fig. 09(c)), 
quantitatively it differs slightly due to the orthogonality 
of the field- and anisotropy directions in this case. 



The model does not directly account for the addition- 
ally observed high field transitions; it predicts continuous 
evolution toward full saturation. However, the model as- 
sumes a distortion-induced finite anisotropy D(B), while 
the associated elastic energy cost is not included in the 
all-magnetic Hamiltonian. The amount of magnetic en- 
ergy the system gains upon having a distortion (and thus 
finite D) can be approximated by taking the difference 
between the energy of an isotropic spin configuration (e.g. 
a canted 120° configuration, with all three spins tilted 
away from the field-direction, while their projections in 
the orthogonal plane keep mutual 120° angles) and the 
minimum energy solution of eq. [2] (lower inset of Fig. [4]). 

The first order transitions at B^ 5 and B^ 3 can then be 
identified as the point where the magnetic energy gain no 
longer outweighs the elastic cost of the distortion, upon 
which the system reverts to the undistorted triangular 
lattice, which is corroborated by the observed isotropy in 
M above b\§ (Fig. [2|). We note the estimated magnetic 
energy gain at this point is ~ 0.42 meV per unit cell (3 
spins), which is approximately the temperature scale of 
the experimental data (3kT at 1.5 K is ~ 0.39 meV). This 

is consistent with the fact that b\§ and B^ are observed 
to shift toward lower fields with increasing temperature. 
For an isotropic lattice above £?J 5 , the model predicts a 
degenerate set of spin structures, among which the afore- 
mentioned canted 120° structure. 

Despite its satisfactory and intuitive results, our simple 
model has its limitations. Though the low field collinear 
phases can be modeled using eq. [TJ being a phenomeno- 
logical model meant to describe the high field phases, it 
does not capture the complex helical ferroelectric phase. 



A full quantitative description of CuFe02 would require 
the inclusion of finite temperature, three dimensional- 
ity (recent work showed the significance of interplane 
coupling a 14 i 15 i 29 i 3Q ), a more refined phonon model, quan- 
tum spins and possibly other interactions^. 

Concluding, through pulsed field magnetization ex- 
periments, the metamagnetic staircase characteristic of 
CuFe02 was extended to up to 58.3 T, revealing an ad- 
ditional first order phase transition for both magnetic 
field configurations, which is proposed to be due to a 
reversed spin Jahn Teller transition. Above this transi- 
tion, virtually complete isotropic behavior is retrieved. 
A highly consistent phenomenological rationalization for 
the magnetization process in both magnetic field con- 
figurations is developed, combining for the first time all 
magnetic terms deemed of importance in CuFe02- Nu- 
merical simulations based on the corresponding classical 
model prove the pertinence of both spin-lattice and field- 
dependent anisotropy interactions in CuFe02- Combined 
with the magnetization measurements, a recovery of the 
undistorted triangular lattice structure is anticipated at 
high fields. The underlying intuitive concept of progres- 
sive symmetry increase as the degree of frustration in 
spin Jahn- Teller distorted systems diminishes is rather 
universal, as it relies solely on energy arguments. In- 
deed, a similar high field transition has been observed 
recently in HgC^CU and corresponding transitions may 
be expected in related spinel systems^. 
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